Exactly solvable model of uniaxial ferroelectrics 
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An exactly solvable lattice model with infinite-range potential is applied to uniaxial ferroelectrics. 
Asymptotically exact expression for free energy as a function of an order parameter at any temper- 
atures is obtained. Effect of thermal expansion of lattice unit cell is taken into account. The free 
energy expansion in powers of order parameter in the vicinity of critical point is presented. Cor- 
rections to Landau expansion are obtained. In particular, it is shown that summand with external 
field contains a contribution of higher powers over order parameter. 
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INTRODUCTION 

Phenomenological description of ferroelectrics in the 
vicinity of critical point outside the Ginzburg critical re- 
gion based, as a rule, on the Landau phenomenological 
theory [Tj. In this theory, the existence some "order pa- 
rameter" ■0 postulated with additional assumption on free 
energy F (ip) expansion in powers of order parameter in 
the vicinity of critical point Tc 
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where a (T) , b (T) , . . . are some phenomenological coef- 
ficients depending on temperature T. Phase transition 
in this theory due to the a (T) sign change in the critical 
point: 



a{T) ^ a{T ~Tc) , a> 0. 
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The other coefficients b (T) , b (T) in sufficiently small 
vicinity Tc assumed as some positive constants. 

The main problems in this approach, beyond the Lan- 
dau expansion justification, is determination of the co- 
efficients a{T), b{T), c{T), .... Theoretical determi- 
nation of these coefficients is closely connected to the 
justiffcation of the Landau expansion problem and es- 
tablishment connection between the system Hamiltonian 
and these coefficients (in some analogy with well known 
Mayer expansion in statistical thermodynamics of non- 
ideal gases) . Determination of the Landau expansion co- 
efficients using experimental data encounters difficulties 
due to narrowness of the applicability range of the ex- 
pansion 
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where Gi is the Ginzburg parameter (size of the critical 
region) . 

The aim of this paper is application of the exactly solv- 
able lattice model with infinite-range interatomic poten- 
tial [2] to uniaxial ferroelectrics. For this model the exact 



expression for free energy as function of temperature and 
Hamiltonian parameters is obtained. Temperature de- 
pendence of dipole moment mean value is deduced. Exact 
expression for free energy as function of order parameter 
at arbitrary temperatures is obtained and exact Landau- 
like expansion in vicinity of the Curie point is derived. It 
is shown that exact expansion for this model contains es- 
sential corrections to the Landau expansion. Comparison 
between theory and experimental data is discussed. 



THE MODEL AND PARTITION FUNCTION 

Uniaxial ferroelectrics has a distinguished direction for 
dipole moments orientation in the system; hence the 
dipole moments projections have two possible values ±P. 
Hamiltonian of this model suppose in following form 
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where is Pj a dichotomous variable with values ±P, asso- 
ciated with j-th site, —J/NPjPk is the spacing indepen- 
dent interactions between electrical dipoles Pj and Pk , E 
is an external electric field strength. Thus, in this model 
dipoles interaction energy does not depend on distance 
between these dipoles. 

Partition function of the system is 
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where /3 = 1/T. 

Let us introduce the order parameter as mean value of 
a dipole moment (P) : 
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Using well known Stratonovich-Hubbard identityf? un- 
der the condition Re a > 
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and appropriate change of variable calculation of the par- 
tition function (|5| reduced to simple quadrature: 

+00 

Z{N,E) = -^ J exp{-Nfiy,E)}dy, (8) 



where 
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Integral Q converges by any values of the parameters 
TV, /3, J, P, E, but unfortunately its analytical calculation 
in elementary functions is impossible. 

To perform analytical evaluation of the partition func- 
tion Q in thermodynamical limit 3> 1, let us expand 
the function ^ in powers of y: 
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To calculate this integral in case (l — 13 JP^) > it is 
enough to keep in expansion ( |10[ ) the terms up to second 
power of y; otherwise, if (l — fSJP'^) < 0, it is neces- 
sary to keep the terms at least up to fourth power of 
y. Thus, the critical temperature is defined by the 
equaUty (l - /3 JP^) = Q: 



Tc = JP. 

Let us introduce dimensionless temperature r 
Then we have 
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where £ is dimensionless external field equal to ratio of 
the dipole energy EP in this external field E to the crit- 
ical temperature Tc- 
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The minimum point of function / (y, £) over variable 
y obey the following equation: 

^:^=r(,-^)-tanh, = 0. (14) 



Note, in critical point r = 1 the second derivative 
fyy {y, £) in presence of an external field is nonzero 
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therefore integral (|8| allows evaluation by means the 
Laplace method. As a result we have in thermodynamic 
limit 
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where yo is solution of the equation (14 1. Note, the vari 



ables yo and £ are interconnected by the equation (14 1. 



ORDER PARAMETER AND FREE ENERGY 

Let us calculate the mean value of the dipole mo- 
ment (P) Q. With account of (14 1, we have: 

^ ^p df{yo,£) 

£\^ '^^ (17) 
rPlyo ] = P tanh yo ■ 
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This formula establishes connection between order pa- 



rameter (P) and point of minimum of function (12 1 
Then free energy per dipole has the following form 
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where yo is solution of equation (14), connected with or 



der parameter by relation (17 1 



Let us express yo via order parameter (P) and substi 

' P 



tute the result into (18 1 
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As a result, we have asymptotically exact expression for 
free energy 74((P)) via order parameter (P): 



A((P),£-)=T 




(20) 



In terms of dimensionless order parameter P^ ~ 
(P) /P (— 1 < cr < 1) the free energy has the following 
form: 



A{cT,£)^T\ 





















[2 




vl-^JJ 



- 5 In (1 + a) 



(21) 



£ In (1 - a) 



3 



(quadratic with respect to external field term omitted be- 
cause it does not depend on order parameter). Hence, for 
systems with Hamiltonian Q the law of corresponding 
states takes place, i.e. equations of state of these systems 
in dimensionless variables t, tr, £ have an identical form. 



LANDAU-LIKE EXPANSION 



Expansion of free energy ( 21 1 in vicinity of the critical 



point has the following form: 
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There are two essential differences between this expan- 
sion and theory Landau. 

1. In contrast to Landau expansion, all the coefficients 



in (22 1 depend on temperature. Moreover, the coef- 



ficient at even powers change their signs with tem- 
perature decreasing. This fact is not essential if 
1 — r ^ 1, but it prevented coefficients determina- 
tion by experimental data beyond small vicinity of 
critical point. Hence, in general out of small vicin- 
ity of critical point it should take into account not 
only the higher powers of a terms, but also change 
of the coefficients with temperature. 



2. The term with external field in (22) contains not 



only a linear over order parameter a contribution, 
but also a power series over a. 



SPONTANEOUS POLARIZATION 

Spontaneous polarization can be found minimum free 
energy (21 1 condition at £ = 0: 
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A schematic graph of function (r) presented on fig. [T] 
To perform a comparison between long-range inter- 
actions uniaxial ferroelectrics model with experimental 
data we should take into account some additional fac- 
tors. 

1. We should return to dimensional polarization Pg 
and save dimensionless temperature r. 

2. We should to take into account thermal expansion 
of materials and related effects. The point is that 
thermal expansion leads to unit cell deformation, 
effective ionic charges change etc. 
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Figure 1. Dependence of dimensionless spontaneous polar- 
ization CT on dimensionless temperature t of uniaxial ferro- 
electrics with long-range potential. 



Note, the spontaneous polarization Pg related to dipole 
moment mean value (5) and specific volume (i.e. volume 
per dipole moment) i7 via following formula 



P. = 



(24) 



Let us find the ratio Pg (r) to (j{t) with account defi- 
nition a = (P) /P: 
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where I oc fi^/"^ is a typical linear size of dipole, q* is a 
typical value of the effective ionic charges in dipoles. 

Change of the temperature leads to change of I (ther- 
mal expansion) and q* (change of the ionic contribution 
of chemical bonds due to redistributions of the electrons) . 
In general, both of values I and q* are increasing functions 
of temperature, bur there is no any certain information 
on the ratio jj" temperature behavior. 

Let us consider the experimental data on tempera- 
ture dependence of spontaneous polarization P^ (r) of the 
uniaxial ferroelectrics triglycine sulfate and theoretical 
curve for dimensionless polarization (t(t) for uniaxial fer- 
roelectrics with long-range interatomic potential. These 
results presented on Fig. [2] 

The temperature dependence of the ratio P, (t) / a (t) 
presented on Fig. [3] 

The ratio Ps (t) /a (r) raising with temperature in- 
crease caused, probably, by the change of ion-covalent 
chemical bond to pure ionic chemical bond and related 
changes of the effective ionic charges in ferroelectrics. 



4 



3,5 

3 ■ 
2,5 ■ 

2 
1,5 ■ 

1 ■ 
0,5 



tain 



-Ps, |jcrm2 



0,84 0,! 



88 0,9 0,92 0,94 0,96 0,' 



Figure 2. Dependence of triglycine-sulfate spontaneous polar- 
ization Ps on dimensionless temperature r j3j (upper curve) 
and theoretical dependence of dimensionless polarization o"(r) 
for long-range interactions model of uniaxial ferroelectrics 
(lower curve). 
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Figure 3. Dependence of the ratio Ps (r) /a (r) for triglycine- 
sulfate. 



Thus, in general case variations of both the unit cell ge- 
ometry and the effective ionic charges due to the materi- 
als thermal expansion should be considered. 



EFFECTS OF THE THERMAL EXPANSION 



Let us express free energy ( 22 1 via polarization Pg , 



defined by (24) 
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where fl (r) and P (r) are the unit cell volume and the 
dipole electric moments as functions of temperature t, 
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Spontaneous polarization P^ in long-range uniaxial fer- 
roelectrics model can be find by this expression mini- 
mization at any values of temperature. In vicinity of the 
critical point |1 — t| ^ 1 the free energy expansion in 



powers of the polarization P, can be obtained from (22 1 
with replacing a on | ^[t\ | As a result, we obtain 
the following coefficients of this expansion: 
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respectively. Substituting this relation into (21 1 we ob- 



CONCLUSION 

The exact solution of long-range uniaxial ferroelectrics 
model with external electric field presence is obtained. 
Effect of thermal expansion is taken into account. The 
compact expression for free energy as function of order 
parameter, temperature, and external electric field is de- 
rived. It is shown that coefficients of the Landau-like 
expansion in critical point vicinity have essential temper- 
ature dependence, up to their signs change. It is shown 
that thermal expansion of unit cell leads to change of the 
elementary dipole moments. The last leads to additional 
variation of the expansion coefficients. The external field 
term in free energy contains not only linear with respect 
to order parameter summand, but some infinite series in 
odd powers of order parameter. 
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